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Nanotube patterning may occur naturally upon the spontaneous self-assembly of biomolecules 
onto the surface of single- walled carbon nanotubes (SWNTs). It results in periodically alternating 
bands of surface properties, ranging from relatively hydrophilic to hydrophobic, along the axis of 
the nanotube. Single Chain Mean Field (SCMF) theory has been used to estimate the free energy 
of systems in which a surface patterned nanotube penetrates a phospholipid bilayer. In contrast 
to un-patterned nanotubes with uniform surface properties, certain patterned nanotubes have been 
identified that display a relatively low and approximately constant system free energy (< ±10 kT) 
as the nanotube traverses through the bilayer. These observations support the hypothesis that 
the spontaneous self-assembly of bio-molecules on the surface of SWNTs may facilitate nanotube 
transduction through cell membranes. 
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Understanding the mechanism of transduction of car- 
bon nanotubes through cell membranes is a challeng- 
ing undertaking from many perspectives. Apart from 
fundamental interest, details of the mechanism may an- 
swer questions about the cytotoxicity of nanotubes 1-4 
and their potential use as delivery vehicles^— or to 
probe bilayer propertie a n i 12 . However, despite consid- 
erable effort devoted to this question, a consensus on the 
mechanism has not yet been reached and direct evidence 
of spontaneous translocation of nanotubes through the 
membranes of cells is still lacking. 

In a recent paper, Pogodin and Baulin 13 considered the 
thermodynamics of a system in which an uncharged nan- 
otube with uniform surface properties penetrates a phos- 
pholipid bilayer. A coarse grained model of the phos- 
pholipid molecule was adopted, which had been previ- 
ously shown to adequately characterise the key thermo- 
dynamic properties of a phospholipid bilayer in a fluid 
phased. Estimates were made of the free energies of 
the equilibrium states for the system as the nanotube 
was progressively moved perpendicularly into the bilayer 
using a numerical implementation of the Single Chain 
Mean Field (SCMF) theory! 4 .. The SCMF methodology 
and three models of the phospholipid bilayer have been 
discussed in details in Ref. .141 The simplest 3-beads 
model of phospholipids has been proved to be successful 
in describing the thermodynamic properties of the bi- 
layer. In essence it involves a coarse grained description 
of the lipid molecule where the monomers are grouped 
into two types of beads, one hydrophilic representing the 
polar heads and two hydrophobic representing the tails 
of the lipid. The size and the interaction parameters of 
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the beads we adjusted to reproduce the essential thermo- 
dynamic properties of a fluid phase of the phospholipid 
bilayer such as the thicknesses of the layer and the hy- 
drophobic core, the equilibrium area per lipid and the 
compressibility constant. 

The objective of these calculations was to determine 
whether nanotubes of various diameters and surface 
properties might penetrate a bilayer as a consequence 
of their thermal motion. Thus, nanotubes of 1, 2.43 and 
4.86 nm diameter were considered, characterised by an 
energy per contact with the coarse grained phospholipid 
tail, Et, ranging from kT (representing steric repulsion) 
to —6.3 kT, which corresponds to strong hydrophobic at- 
traction. The perpendicular orientation was chosen since 
this represents the minimum contact area between nan- 
otube and phospholipids per unit depth of penetration 
and hence the minimum free energy of interaction. An 
output of the SCMF calculations was the equilibrium free 
energies and the spatial mean field concentration distri- 
bution of the phospholipid heads and tails in the bilayer, 
which varied as a consequence of nanotube penetration. 
The model thus demonstrated the structural rearrange- 
ment of phospholipids at the molecular level that was 
induced by insertion of the nanotube, and the equilib- 
rium free energy change of the system for each position 
of the nanotube. 

In summary, the calculations showed that the free en- 
ergy change of the systems for et = —2.1 kT rose mono- 
tonically with increasing nanotube penetration and were 
substantial at full penetration (e.g. for the 2.43 nm di- 
ameter nanotube the free energy of the system at full 
penetration was about 100 kT). For et = —4.2 kT, and 
particularly for et = —6.3 kT, the initial penetration of 
the nanotube resulted in a significant fall in the system 
free energy (e.g. to ca. —80 kT for penetration of a 2.43 
nm diameter nanotube with et — —6.3 kT to the centre 
of the bilayer) . Further penetration of these nanotubes 
then led to a steep rise in the system free energy. Overall, 
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FIG. 1. (A) Possible patterns of carbon nanotubes. The dense colours indicate the magnitude of the interaction parameter 
with the phospholipid tails, et, for the corresponding section of nanotube. The translucent colouring is provided to indicate the 
interaction range where the phospholipids feel the attraction. (B) Concentration profiles of the phospholipid bilayer interacting 
with differently patterned nanotubes. Complete set of concentration profiles is available in supplementary materials section. 



the calculations showed that hydrophilic and weakly hy- 
drophobic nanotubes face a substantial energy barrier to 
penetration, whilst intermediate and strongly hydropho- 
bic nanotubes penetrate little or become entrapped in a 
free energy well within the bilayer. 

Inferences can be drawn from these calculations for the 
transduction of cylindrical nano-objects such as single- 
walled carbon nanotubes (SWNTs) through the outer 
membrane of cells. Untreated SWNTs are significantly 
hydrophobic and their thermal motion might lead to their 
accumulation within the core of the cell membrane, but 



they are unlikely to translocate across the membrane in 
view of the steep energy barrier that they face to pass 
out of the membrane core. Hydrophilically functionalised 
SWNTs, for example pegylated SWNTs^ face a substan- 
tial energy barrier in even penetrating the membrane to 
its core. These thermodynamic calculations create a co- 
nundrum since numerous experimental reports exist that 
show the accumulation of SWNTs within the cytoplasm 
of cell a 16 ! 17 . How may this be? 

Unlike the model system that is envisaged in Ref. 
Il3l . the environment around a cell is complex and com- 
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Nanotube position (nm) 

FIG. 2. Free energy cost AF versus nanotube position of pat- 
terned SWNTs with diameter 1 nm SI and S2 in comparison 
with uniform nanotubes with different interaction parameters 
with the hydrophobic core of the phospholipid bilayer, et^. 
The unperturbed phospolipid bilayer location is indicated by 
the translucent colouring. 



prises of diverse biomolecular species that might interact 
with SWNTs. Indeed, a number of publications have re- 
ported the ordered self-assembly of polar lipid a 18 i 19 , sin- 
gle stranded DNA2£~— , polysaccharide s 23 ' 24 , amphiphilic 
proteins 2 ^— and even vitamins— onto nanotubes. The 
electrostatics also give rise to the most general patterns 
on the nanotubes^. This self-assembled patterning oc- 
curs spontaneously upon mixing the nanotubes with the 
patterning agent in aqueous solution. The resulting 
molecular structures commonly take the form of discrete 
hydrophilic rings along the axis of the nanotube in the 
case of polar lipids, to helices for polysaccharides and 
DNA. Thus, in practice, in cell culture systems nanotubes 
may not have homogeneous surface properties but may 
display a distinct regular patterning. Furthermore, it is 
not evident that a nanotube that has been naturally pat- 
terned in this way interacts with phospholipid bilayers 
in a similar manner to a naked nanotube or to a nan- 
otube with a homogeneous adsorption layer. Our point 
of view is supported by experimental evidence that or- 
dered arrangements of hydrophilic and hydrophobic sur- 
face functional groups can alter the penetration of spher- 
ical nanoparticles through the cell membranes 3 ^. 

Consider then a patterned hydrophobic nanotube. 
Along the nanotube the self-assembly of polar lipids (or 
other biomolecular species) leads to equi-spaced rings of 
different relative hydrophobicity to the naked nanotube. 
Moving along the axis of the nanotube, the surface char- 
acteristics alternate periodically as indicated in Figure 
QJA). The spatially segregated surface characteristics of 
such a nanotube is thus substantially different from the 
uniform surface character assumed in the previous SCMF 
calculations—, and this may influence the free energy 
change of the system when the nanotube penetrates a 



phospholipid bilayer. Note, that the surface patterning 
has bigger effect on the translocation through a bilayer 
than the shape or geometry of nanoparticles^i, especially 
for small particles. 

Comparison of the energy curves of uniform nanotubes 
at different positions^ suggest that alternation of stripes 
with certain interaction can, in principle, reduce consid- 
erably the energy barrier of translocation. For example, 
two stripes with opposite energies placed together may 
cancel the contribution of each other. Since the thick- 
ness of the hydrophobic core is about 2 nm, the width of 
alternating stripes in this case should be of order 1 nm 
so that the core is in contact simultaneously with two 
opposite stripes. Two distinct patterns were considered 
for a 1 nm diameter nanotube (Figure [HA) , SI and S2) 
that differ only in the relative widths of the two sets of 
rings (i.e. whereas B=1.00 nm for SI, D=0.90 nm and 
E=1.10 nm for S2). Both SI and S2 are characterised 
by alternating rings of interaction energy St = —2.1 and 
—6.3 kT respectively. A larger nanotube of diameter 2.43 
nm was also considered. In this case five different pat- 
terns were considered (Figure QJA) , LI, L2, L3, L4 and 
L5). The effect of the bottom of the nanotube was also 
investigated, L4 and L5 only differ in the interaction of 
the bottom part, while PI has an edge with the same 
interaction parameter as L4. The corresponding concen- 
tration profiles of nanotubes with similar patterns, S2, 
PI, L4 and L5, are shown in Figure [IjB). We emphasise 
that these patterns are purely representative and their 
dimensions and energies are not those for any specific 
biomolecular system. 

The concentration profiles for the patterns LI, L2, L3 
look identical, while the energy of penetration of these 
patterns is quite different. This is because the same po- 
sitions of the lipids around the nanotube may lead to 
different enthalpic contributions. Thus, the snapshots 
usually provided by MD simulations may not be suffi- 
cient to distinguish between different scenarios, while the 
equilibrium energy of insertion may be crucial for under- 
standing the mechanisms of insertion into phospholipid 
bilayers. 

The free energy of penetration of these variously pat- 
terned nanotubes into a coarse grained representation of 
a phospholipid bilayer was estimated using a numerical 
implementation of SCMF, in the manner described by 
Pogodin and Baulini 3 .. The resulting free energies for 
the small nanotubes are shown in Figure [2j For the nan- 
otube pattern SI with equi-sized rings (Figure [T{ A), SI) 
the system free energy always lies within the range to 
—20 kT. Neglecting the end of the nanotube, the mean 
interaction energy along the length of SI is —4.2 kT. 
Comparing SI with the free energy for penetration of 
a nanotube with a uniform surface of interaction energy 
E4, -4.2 kT (Figure^ E4) it is evident that no steep en- 
ergy barrier is encountered by the patterned nanotube. 
Variation of the relative sizes of the rings (to give Figure 
QJA), S2) results in the free energy of the system being 
always within the range 0< ±10 kT, allowing relatively 
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FIG. 3. Free energy cost AF versus nanotube position of SWNTs with diameters 2.43 nm with different pattern, LI, L2, 
L3 (A) and different end face L4, L5, PI (B) in comparison with uniform nanotubes with different interaction parameters 
with the hydrophobic core of the phospholipid bilayer, £t-£. The unperturbed phospolipid bilayer location is indicated by the 
translucent colouring. 



TABLE I. Energy cost of full insertion, AFf u u, the maximal amplitude (positive or negative) of the insertion free energy, 
AF max , insertion distance corresponding to the minimum energy, d m i n , maximal force and maximal pressure, for piercing the 
DMPC phospholipid bilayer. 
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unimpeded passage of the nanotube through the bilayer. 
Note, that the oscillations of the energy curve is due to 
abrupt passage from one stripe to another. However, one 
can expect that a helical pattern would allow for smooth 
transition from one minimum to another similar to screw- 
driving. Since the error bar of these curves is about few 
kT, this pattern implies a zero energy cost of transloca- 
tion. 

A similar behaviour was observed for the larger diam- 
eter nanotube (Figure [^A), LI, L2 and L3). Judicious 
choice of patterning (L3) resulted in free energies for pen- 
etration within the range —20 to +10 kT with no steep 
energy barrier. Injudicious patterning (LI and L2) re- 
sulted in free energy profiles that varied significantly with 
depth of penetration and which both presented very sig- 
nificant energy barriers to penetration. The result for 
pattern LI suggests a critical dependence upon the in- 
teraction energy of the end face of the nanotube since 
the pattern LI differs from that of a smaller nanotube 
SI only insofar as the end face has been changed to the 
highest hydrophobicity (c.f. —6.3 kT in LI and —2.1 
kT in SI). The free energy profile for LI is thus similar 
in characteristic to the uniformly hydrophobic nanotube, 
though presents a significantly steeper energy barrier to 
full penetration. For L2 the effect is opposite, the lower 



hydrophobicity of the end face and rings causes the free 
energy profile to resemble that for a nanotube of uni- 
formly low hydrophobicity and again presents a steep 
energy barrier at a relatively low extent of penetration. 
Finally, the error bars on the point at nm for curve L3 
in Figure A) are included to show the standard error 
of six separate SCMF calculations for this arrangement 
and extent of penetration (approximately 10 kT). The 
characteristic values of the energies, forces and pressures 
for different patterns are summarized in Table HI 

The effect of the end face was investigated in Figure 
(3{B). The pattern of L4, L5 and PI is similar to a suc- 
cessful patterning of the small nanotube S2, only the end 
face is different. The larger nanotube L5 has the same 
end face as S2 and shows no steep energy barrier for pene- 
tration, although the energy is slightly shifted to positive 
values. The same pattern, but different end face, L4, re- 
sults in serious changes in the penetration energy in the 
beginning of insertion. In fact, the energy curve for L4 
follows the curve for homogeneous nanotube E6 which 
has the same interaction parameter as the end face. In 
turn, when the nanotubes are fully inserted, the end face 
does not influence and the curves L4 and L5 coincide. 
Modification of the shape of the end of the nanotube 
does not lead to serious changes. Since our calculations 
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consider equilibrium insertion, only the area of contact 
of different stripes matters and the nanotube with sharp 
tip PI behaves similar to L4. However, one can expect 
that the shape of the tip may be important for dynamics 
of the insertion. 

In the calculations reported here, no attempt has been 
made to determine the optimum surface patterning along 
the axes of the nanotubes that results in the most uniform 
system free energy throughout bilayer penetration that 
we leave for future work. Nor have we sought to deter- 
mine the optimal sets of et's. Rather, calculations have 
been conducted for discrete patterns of combination of 
previously reported ex's for homogeneous nanotubes^. 
Nevertheless, for the systems studied the effect of nan- 
otube patterning had a significant effect upon the system 
free energy change during bilayer penetration and some 
patterns were shown to result in a relatively uniform free 



energy throughout penetration. In the most preferable 
case, S2, any free energy barrier to penetration is < 10 
kT, which is comparable in magnitude to the standard 
error in the SCMF free energy calculation. Extrapolat- 
ing these observations to the practical case of cell mem- 
brane penetration by SWNTs, it is tempting to specu- 
late that patterning of the tubes by one or other of the 
biomolecules that are commonly present in cell culture 
supernatants may significantly enhance the possibility of 
their transduction through cell membranes. 
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